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Abstract 

The structural, electronic and magnetic properties of small Fe m Rh n clusters having N = m+n < 
8 atoms are studied in the framework of a generalized-gradient approximation to density-functional 
theory. For JV = m|n < 6 a thorough sampling of all cluster topologies has been performed, 
while for N = 7 and 8 only a few representative topologies are considered. In all cases the en- 
tire concentration range is systematically investigated. All the clusters show ferromagnetic-like 
order in the optimized structures. As a result, the average magnetic moment per atom JJ N in- 
creases monotonously, which is almost linear over a wide range of concentration with Fe content. 
A remarkable enhancement of the local Fe moments beyond 3 \xb is observed as result of Rh dop- 
ing. The composition dependence of the binding energy, average magnetic moment and electronic 
structure are discussed. 

PACS numbers: 75.75.+a, 36.40. Cg, 75.50.Bb, 73.22.-f 
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I. INTRODUCTION 



Alloying elements with complementary qualities in order to tailor their physical behavior 
for specific technological purposes has been a major route in material development since the 
antiquity. Cluster research is no exception to this trend. After decades of systematic studies 
of the size and structural dependence of the most wide variety of properties of monoelement 
particles, the interest has actually been moving progressively over the past years towards 
investigations on finite-size binary alloys.- The magnetism of transition-metal (TM) clusters 
opens numerous possibilities and challenges in this context.- - — For example, one would 
like to understand how to modify the magnetic characteristics of clusters, in particular the 
saturation magnetization and the magnetic anisotropy energy (MAE), as it has been done 
in solids. This would indeed allow one to design new nanostructured materials from a 
microscopic perspective. Nevertheless, it also true that controlling composition, system size, 
and magnetic behavior sets serious difficulties for both experiment and theory. 

Pure TM clusters such as Fe^y, Co at and Ni^r show spin moments, orbital moments, and 
MAEs that are enhanced with respect to the corresponding periodic solids.— ~— Still, the 
possibilities of optimizing the cluster magnetic behavior by simply tuning the system size 
have been rather disappointing, particularly concerning the MAE, which remains relatively 
small — despite being orders of magnitude larger than in solids^- — due to the rather weak 
spin-orbit (SO) coupling in the 3d atoms. This is one of the motivations for alloying 3d 
TMs with Ad and 5d elements which, being heavier, are subject to stronger SO interactions. 
In this context it is useful to recall that large nanoparticles and three dimensional solids of 
these elements are non-magnetic. However, at very small sizes the 4d and 5d clusters often 
develop a finite spontaneous low-temperature magnetization, due to the reduction of local 
coordination and the resulting <i-band narrowing.— ~— The first experimental observation of 
this important finite-size effect has been made by Cox et al. by performing Stern-Gerlach- 
deflection measurements on Rh^v clusters. In this work the average magnetic moments per 
atom ~p N = 0.15-0.80yU£ have been experimentally determined for N < 30-50 atoms.— In 
view of these contrasting features one expects that 3d-4d and 3d-5d alloy clusters should 
show very interesting structural, electronic and magnetic behaviors. 

The purpose of this paper is to investigate the ground-state properties of the small FeRh 
clusters in the framework of Hohenberg-Kohn-Sham's density functional theory— Besides 
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the general interest of the problem from the perspective of 3d-4d nanomagnetism, these 
clusters are particularly appealing because of the remarkable phase diagram of FeRh bulk 
alloys.— In the case of Fe5oRh 50 the magnetic order at normal pressure and low temperatures 
is antiferromagnetic (AF). As the temperature increases this a" phase undergoes a first order 
transition to a ferromagnetic (FM) state, the a' phase, which is accompanied by a change 
in lattice parameter. The corresponding transition temperature Tf a " increases rapidly 
with increasing external pressure P, eventually displacing the FM a' phase completely for 
P > 7 GPa (Tf a " ~ 290 K for FesoRhso at normal pressure). Moreover, T"' a " decreases 
very rapidly with decreasing Rh content. At low pressures the FM a' phase undergoes a FM 
to paramagnetic (PM) transition at (Tq — 670 K).— In addition, the properties of a- FeRh 
bulk alloys have been the subject of first principles and model theoretical investigations.— 
In particular these show that the relative stability of the FM and AF solutions depends 
strongly on the interatomic distances. Such remarkable condensed-matter effects enhance 
the appeal of small FeRh particles as specific example of 3d-4d nanoscale alloy. Investigations 
of their magnetic properties as a function of size, composition, and structure are therefore 
of fundamental importance. 

The remainder of the paper is organized as follows. In Sec. |TT] the main details of the 
theoretical background and computational procedure are presented. This includes in par- 
ticular a description of the strategy used for exploring the cluster energy landscape as a 
function of geometrical conformation and chemical order. The results of our calculations for 
FeRh clusters having N < 8 atoms are reported in Sec. II I II by analyzing the concentration 
dependence of the cohesive energy, the local and average magnetic moments, and the spin- 
polarized electronic structure. Finally, we conclude in Sec. IIVI with a summary of the main 
trends and an outlook to future extensions. 

II. COMPUTATIONAL ASPECTS 

The calculations reported in this work are based on density functional theory— as imple- 
mented in the Vienna ab initio simulation package (VASP).— The exchange and correlation 
energy is described by using both the spin-polarized local density approximation (LDA) 
and Perdew and Wang's generalized-gradient approximation (GGA).— The VASP solves 
the spin-polarized Kohn-Sham equations in an augmented plane-wave basis set, taking into 
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account the core electrons within the projector augmented wave (PAW) method.— This is 
an efficient frozen-core all-electron approach which allows to incorporate the proper nodes 
of the Kohn-Sham orbitals in the core region and the resulting effects on the electronic 
structure, total energy and interatomic forces. The 4s and 3d orbitals of Fe, and the 5s 
and 4d orbitals of Rh are treated as valence states. The wave functions are expanded in a 
plane wave basis set with the kinetic energy cut-off E max = 268 eV. In order to improve the 
convergence of the solution of the selfconsistent KS equations the discrete energy levels are 
broadened by using a Gaussian smearing a = 0.02 eV. The validity of the present choice 
of computational parameters has been verified.— The PAW sphere radii for Fe and Rh are 
1.302 A and 1.402 A, respectively. A simple cubic supercell is considered with the usual 
periodic boundary conditions. The linear size of the cell is a = 10-22 A, so that any pair 
of images of the clusters are well separated and the interaction between them is negligible. 
Since we are interested in finite systems, the reciprocal space summations are restricted to 
the T point. 

Although the potential advantages of alloying magnetic 3c? elements with highly-polarizable 
4d or 5d elements are easy to understand, the problem involves a number of serious practical 
challenges. Different growth or synthesis conditions can lead to different chemical orders, 
which can be governed not just by energetic reasons but by kinetic processes as well. For 
instance, one may have to deal with segregated clusters having either a 4d core and a 3d 
outer shell or vice versa. Post-synthesis manipulations can induce different degrees of in- 
termixing, including for example surface diffusion or disordered alloys. Moreover, the inter 
atomic distances are also expected to depend strongly on size and composition. Typical 
TM-cluster bond-lengths are in fact 10-20% smaller than in the corresponding bulk crystals. 
Taking into account that itinerant 3d-electron magnetism is most sensitive to the local and 
chemical environments of the atoms,— >2tM it is clear that controlling the distribution of the 
elements within the cluster is crucial for understanding magnetic nanoalloys. 

Systematic theoretical studies of binary-metal clusters are hindered by the diversity of 
geometrical conformations, ordered and disorder arrangements, as well as segregation ten- 
dencies that have to be taken into account. This poses a serious challenge to both first- 
principles and model approaches. In order to determine the interplay between cluster struc- 
ture, chemical order and magnetism in FeRh clusters we have performed a comprehensive 
set of electronic calculations for clusters having N < 8 atoms. In the present paper we focus 
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on the most stable cluster structure and magnetic configuration, which are determined by 
exploring the ground-state energy landscape.— This is a formidable task, since one needs 
to consider a large, most possibly complete and unbiased set of initial structures. Such a 
thorough geometry optimization must include not only the representative cluster geometries 
or topologies, but also all relevant chemical orders. This requires taking into account all 
distributions of the Fe and Rh atoms for any given size and composition. These two aspects 
of the problem of determining the structure of nanoalloys are discussed in more detail in the 
following. 

The different cluster topologies are sampled by generating all possible graphs for N < 6 
atoms as described in Ref. 



38 (see also Ref. 



4 if ). For each graph or adjacency matrix it is 



important to verify that it can be represented by a true structure in D < 3 dimensions. 
A graph is acceptable as a cluster structure, only if a set of atomic coordinates Ri with 
% — 1, . . . , N exists, such that the interatomic distances R^ satisfy the conditions Rij = Rq 
if the sites i and j are connected in the graph (i.e., if the adjacency matrix element = 1) 
and R^ > Ro otherwise (i.e., if Ay = 0). Here Rq refers to the nearest neighbor (NN) 
distance, which at this stage can be regarded as the unit of length, assuming for simplicity 
that it is the same for all clusters. Notice that for N < 4 all graphs are possible cluster 
structures. For example, for N = 4, the different structures are the tetrahedron, rhombus, 
square, star, triangular racket and linear chain.— However, for N > 5 there are graphs, 
i.e., topologies, which cannot be realized in practice. For instance, it is not possible to 
have five atoms being NNs from each other in a three dimensional space. Consequently, for 
N > 5 there are less real structures than mathematical graphs. The total number of graphs 
(structures) is 21 (20), 112 (104), and 853 (647) for N = 5, 6, and 7, respectively.-^ 

For clusters having N < 6 atoms all these topologies have indeed been taken as starting 
points of our structural relaxations. Out of this large number of different initial configura- 
tions the unconstrained relaxations using VASP lead to only a few geometries, which can 
be regarded as stable or metastable isomers. For larger clusters (N = 7 and 8) we do not 
aim at performing a full global optimization. Our purpose here is to explore the interplay 
between magnetism and chemical order as a function of composition for a few topologies that 
are representative of open and close-packed structures. Taking into account our results for 
smaller sizes, and the available information on the structure of pure Fe^ and Rli^r clusters, 
we have restricted the set of starting topologies for the unconstrained relaxation of FeRh 



heptamers and octamers to the following: bicapped trigonal bipyramid, capped octahedra, 
and pentagonal bipyramid for N = 7, and tricapped trigonal bipyramid, bicapped octahe- 
dra, capped pentagonal bipyramid and cube for N = 8. Although, the choice of topologies 
for N = 7 and 8 is quite restricted, it includes compact as well as more open structures. 
Therefore, it is expected to shed light on the dependence of the magnetic properties on the 
chemical order and composition. 

The dependence on concentration is investigated systematically for each topology of 
Fe m Rh n by varying m and for each size N = m + n < 8, including the pure Fe^ and 
Rh^r limits. Moreover, we take into account all possible non-equivalent distributions of the 
m Fe and n Rh atoms within the cluster. In this way, any a priori assumption on the 
chemical order is avoided. Obviously, such an exhaustive combinatorial search increasingly 
complicates the computational task as we increase the cluster size, and as we move away 
from pure clusters towards alloys with equal concentrations. Finally, in order to perform the 
actual density-functional calculations we set for simplicity all NN distances in the starting 
cluster geometry equal to the Fe bulk valued i?o = 2.48 A. Subsequently, a fully uncon- 
strained geometry optimization is performed from first principles by using the VASP.— The 
atomic positions are fully relaxed by means of conjugate gradient or quasi-Newtonian meth- 
ods, without imposing any symmetry constraints, until all the force components are smaller 
than the threshold 5 me V/A. The convergence criteria are set to 10 5 eV/A for the energy 
gradient, and 5 x 10 -4 A for the atomic displacements.— The same procedure applies to 
all considered clusters regardless of composition, chemical order, or total magnetic moment. 
Notice that the diversity of geometrical structures and atomic arrangements often yields 
many local minima on the ground-state energy surface, which complicates significantly the 
location of the lowest-energy configuration. 

Lattice structure and magnetic behavior are intimately related in TMs, particularly in 
weak ferromagnets such as Fe and its alloys.— On the one side, the optimum structure and 
chemical order depend on the actual magnetic state of the cluster as given by the average 
magnetic moment per atom JJ N and the magnetic order. On the other side, the magnetic 
behavior is known to be different for different structures and concentrations. Therefore, in 
order to rigorously determine the ground-state magnetic properties of FeRh clusters, we have 
varied systematically the value of the total spin polarization of the cluster S z by performing 
fixed spin-moment (FSM) calculations in the whole physically relevant range. Let us recall 
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that S z = (z/f — v±)/2 where v^(y£) represents the number of electrons in the majority 
(minority) states. In practice we start from the non-magnetic state (S™ m = 0) and increase 
S z until the local spin moments are fully saturated, i.e., until the Fe moments in the PAW 
sphere reach /i Fc ~ 4/i B and the Rh moments /iRh — 2.5/i# (typically, S™ ax > 3N/2). The 
above described global geometry optimizations are performed independently for each value 
of S z . These FSM study provides a wealth of information on the isomerization energies, the 
spin-excitation energies, and their interplay. These are particularly interesting for a subtle 
magnetic alloy such as FeRh, and would therefore deserve to be analyzed in some more detail. 
In the present paper we shall focus on the ground-state properties by determining for each 
considered Fe m Rh n the most stable structural and magnetic configuration corresponding to 
energy minimum as a function of S z and of the atomic positions.— 

Once the optimization with respect to structural and magnetic degrees of freedom is 
achieved, we derive the binding energy per atom Eb = [mE(Fe) +nE(Rh) — E(Fe m Rh n )]/N 
in the usual way by referring the total energy E to the corresponding energy of m Fe and 
n Rh isolated atoms. Moreover, for each stationary point of the total energy surface (i.e., 
for each relaxed structure having a nearly vanishing ||Vi?||) we determine the vibrational 




frequencies from the diagonalization of the dynamical matrix. The latter is calculated from 
finite differences of the analytic gradients of the total energy. In this way we can rule out 
saddle points to which the local optimization procedure happens to converge on some occa- 
sions. Only configurations which correspond to true minima are discussed in the following. 
Finally, a number of electronic and magnetic properties — for example, the binding energy, 
the local magnetic moments /ij integrated within the Wigner-Seitz (WS) or Bader atomic 
cells of atom z,— ^ and the spin polarized density of electronic states (DOS) p a (£) — are 
derived from the self-consistent spin-polarized density and Kohn-Sham spectrum. 



III. RESULTS AND DISCUSSION 

In the following we present and discuss results for the binding energy, average and local 
spin moments, and electronic densities of states for N = m + n < 8. 
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FIG. 1. (Color online) Binding energy per atom Eb of Fe m Rh n clusters as a function of the number 
of Fe atoms. The lines connecting the points for each N = m + n are a guide to the eye. 

A. Binding energy and magnetic moments 

In Fig. [1] the binding energy per atom Eb is given as a function of the number of Fe 
atoms m. Besides the expected monotonic increase of Eb with increasing N, an interesting 
concentration dependence is observed. For very small sizes (N < 4) Eb is maximal for 
m = 1 or 2, despite the fact that Eb is always larger for pure Rh than pure Fe clusters. 
This indicates that in these cases the bonding resulting from FeRh pairs is stronger than 
RhRh bonds. Only for m > N — 1, when the number of weaker FeFe bonds dominates, 
one observes that Eb decreases with increasing m. For larger sizes (N > 5) the strength of 
RhRh and FeRh bonds becomes very similar, so that the maximum in Eb is replaced by a 
range of Fe concentrations x = m/N < 0.5 where Eb depends very weakly on m. 

In Fig. |2] the average magnetic moments Ji N of Fe m Rh n are shown as a function of m 
for N < 8. First of all, one observes that ~p N increases monotonously, with the number 
of Fe atoms. This is an expected consequence of the larger Fe local moments and the 
underlying FM-like magnetic order. The average slope of the curves tends to increase with 
decreasing N, since the change in concentration per Fe substitution is more important the 
smaller the size is. The typical increase in ~p, N per Fe substitution is about (1/N)/j,b per 
Fe substitution. Notice, moreover, the enhancement of the magnetic moments of the pure 
clusters in particular for Fe^ (m — N), which go well beyond 3fis, the value corresponding 
to a saturated <i-band in the d 7 s l configuration. In contrast, the moments of pure Rhjv are 
far from saturated except for N = 2 and 7 (see Fig. [2] for m — 0). In this context it is 
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FIG. 2. (Color online) Total magnetic moment per atom JI N of Fe m Rh n clusters as a function of 
number of Fe atoms. The symbols corresponding to each size are the same as in Fig. [TJ The lines 
connecting the points for each N = m + n are a guide to the eye. 
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FIG. 3. (Color online) Local magnetic moment fi a at the Fe and Rh atoms as a function of the 
number Fe atoms m. 

important to recall that a thorough global optimization, for example, by considering a large 
number of initial topologies, could affect the quantitative values of the magnetic moments 
for N = 7 and 8. 

The local magnetic moments in the PAW sphere of the Fe and Rh atoms provide further 
insight on the interplay between 3d and Ad magnetism in Fe m Rh n . In Fig. [3] fip c and /^Rh 
are shown as a function of m for N = 6-8. The Fe moments are essentially given by the 
saturated <i-orbital contribution. For pure Fe clusters the actual values of /xp e within the 
PAW sphere are somewhat lower than 3/i# due to a partial spill-off of the spin-polarized 
density. Notice that the Fe moments increase as we replace Fe by Rh atoms showing some 
weak oscillations as a function of m. The increase is rather weak for a single Rh impurity 
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in FeAr_iRh but becomes stronger reaching a more or less constant value as soon as the 
cluster contains 2 or more Rh (m < N — 2, see Fig. [3]). This effect can be traced back to a d 
electron charge transfer from Fe to Rh which, together with the extremely low coordination 
number, yields a full polarization of the larger number of Fe d holes. On the other side 
the Rh moments are not saturated and therefore are more sensitive to size, structure and 
composition. The values of //Rh are in the range of 1-1.5/xb showing some oscillations as a 
function of m. No systematic enhancement of /iRh with increasing Fe content is observed. 
This behavior could be related to charge transfers effects leading to changes in the number 
of Rh d electrons as a function of m. 

Finally, it is interesting to analyze the role played by magnetism in defining the cluster 
structure by comparing magnetic and non-magnetic calculations. For the smallest FeRh 
clusters (N = 3 and 4) the magnetic energy AE m = E(S Z = 0) — E(Sz) gained upon 
magnetization is higher in the first excited isomer than in the most stable structure. This 
implies that the contribution of magnetism to the structural stability is not crucial, since 
the non-magnetic calculations yield the same ordering, at least concerning the two best 
structures. This suggests that for the smallest sizes the kinetic or bonding energy dominates 
the structural stability, which also explains that the two most stable isomers have different 
topologies. The situation changes for large clusters. For N > 5 one finds a number of FeRh 
clusters for which the optimal structure is actually stabilized by magnetism. For example, 
in Fe4Rh, Fe3Rh 2 , and FeRh 4 the energy ordering of the two most stable isomers would be 
reversed if magnetism were neglected. It should be noted that in these cases the structures 
differ only in the chemical order, not in the topology which is a TBP. In the FeRh hexamers 
the energy differences between the low-lying isomers are more important and only in one 
case, Fe4Rli2, magnetism appears to be crucial for stabilizing the actual optimal structure. 
A similar strong interplay between structure, chemical order and magnetism is expected for 
larger FeRh clusters. 

B. Electronic structure 

It is very interesting to analyze, at least for some representative examples, how the 
electronic structure depends on the composition of magnetic nanoalloys. To this aim we 
report in Fig. H]the spin-polarized ^-electron density of states (DOS) of representative FeRh 
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FIG. 4. (Color online) Electronic density of states (DOS) of FeRh octamers. Results are given 
for the total (solid), the Fe-projected (dotted), and the Rh-projected (dashed) d-electron DOS. 
Positive (negative) values correspond to majority (minority) spin. A Lorentzian width A = 0.02 eV 
has been used to broaden the discrete energy levels. 

octamers having the most stable relaxed configuration among the considered topologies (see 
Sec. [Ill)- Results for pure Feg and Rh 8 are also shown for the sake of comparison. In all the 
clusters, the dominant peaks in the relevant energy range near Bp correspond either to the 
Fe-3d or to the Kh-Ad states. The valence spectrum is largely dominated by these (^-electron 
contributions. In fact the total DOS and the (^-projected DOS are difficult to tell apart. 

First of all, let us consider the DOS of the pure clusters. Our results for Rhg with a 
cube structure are similar to those of previous studies.— They show the dominant <i-electron 
contribution near Ep, with the characteristic ferromagnetic exchange splitting between the 
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minority and majority spin states. In Fig. H] we also included the DOS for Rh 8 with a 
BCO structure, since it allows us to illustrate the differences in the electronic structure 
of compact and open geometries. Moreover, the DOS of pure Rh 8 with BCO structure 
is very useful in order to demonstrate the dependence of DOS on Fe content, since the 
structures of Fe m Rh 8 _ m with m > 1 are similar to the BCO. Both Fe 8 and Rh 8 show 
relatively narrow <i-bands which dominate the single-particle energy spectrum in the range 
— 5eV < e — ef < 3eV. The spin polarization of the DOS clearly reflects the ferromagnetic 
order in the cluster. Putting aside the exchange splitting, the peak structure in the up and 
down DOS p a (e) are comparable. There are even qualitative similarities between the two 
elements. However, looking in more detail, one observes that the effective <i-band width in 
Fe 8 (about 4 eV) is somewhat smaller than in Rh 8 (about 5 eV). Moreover, in Rh 8 the DOS 
at Ef is non- vanishing for both spin directions and the finite-size gaps are very small (see 
Fig. H]). In contrast, the majority <i-DOS is fully occupied in Fe 8 , with the highest majority 
state lying about 0.5 eV below ef- In addition there is an appreciable gap (about 0.1 eV) in 
the corresponding minority spectrum. These qualitative differences are of course consistent 
with the fact that Fe 8 is a strong ferromagnet with saturated moments, while Rh 8 should 
be regarded clS 8b weak unsaturated ferromagnet. 

The trends as a function of concentration reflect the crossover between the previous 
contrasting behaviors. For low Fe concentration (e.g., Fe2Rhg) we still find states with both 
spin directions close to Ef- The magnetic moments are not saturated, although the Fermi 
energy tends to approach the top of the majority band. Moreover, the majority-spin states 
close to Ef have dominantly Rh character. Small Fe doping does not reduce the d-band 
width significantly. Notice the rather important change in the shape of the DOS in Fe2Rh 6 
as compared to the DOS in Rh 8 . This is a consequence of the change in topology from cubic 
to bicapped octahedron (BCO). 

For equal concentrations (Fe4Rh 4 ) the first signs of <i-band narrowing and enhanced 
exchange splitting start to become apparent. The spin-up states (majority band) which in 
Fe2Rli6 contribute to the DOS at ef now move to lower energies (0.3 eV below ef) so that 
the majority band is saturated. Only spin-down (minority) states are found around Ep, 
although there is a significant gap in pi(s) (see Fig. Hj). In the majority band Rh dominates 
over Fe at the higher energies (closer to Ep), while Fe dominates in the bottom of the band. 
In the minority band the participation of Rh (Fe) is stronger (weaker) below Ep and weaker 
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(stronger) above Bp. This is consistent with the fact that the Rh local moments are smaller 
than the Fe moments. 

Finally, in the Fe rich limit (e.g., FeeRh 2 ), the majority-band width becomes as narrow as 
in Fes, while the minority band is still comparable to Rhg. The exchange splitting is large, 
the majority band saturated and only minority states are found close to Bp. As in Fes, Pi{ £ ) 
shows a clear gap at Bp (see Fig. @|. However, the Rh contribution to the minority states 
below ep remains above average despite the relative small Rh content. The Fe contribution 
largely dominates the unoccupied minority-spin DOS, in agreement with the larger local Fe 
moments. 

IV. SUMMARY AND OUTLOOK 

The structural, electronic and magnetic properties of small Fe m Rh n clusters having N = 
m + n < 8 atoms have been investigated systematically in the framework of a generalized 
gradient approximation to density-functional theory. For very small sizes (N < 4 atoms) the 
binding energy Eb shows a non-monotonous dependence on concentration, which implies 
that the FeRh bonds are stronger than the homogeneous ones. However, for larger sizes 
the FeRh and RhRh bond strengths become comparable, so that Eb depends weakly on 
concentration for high Rh content. 

The magnetic order of the clusters having the most stable structures is found to be FM- 
like. Moreover, the average magnetic moment per atom 7% increases monotonously, which 
is almost linear over a wide range of concentration with Fe content. Consequently, the 
energy gain AE m associated to magnetism also increases with the number of Fe atoms. The 
largest part of the spin polarization (about 90%) can be traced back to the local d magnetic 
moments within the PAW sphere of the atoms. The s and p spin polarizations are almost 
negligible in general. A remarkable enhancement of the local Fe moments is observed as 
result of Rh doping. This is a consequence of the increase in the number of Fe d holes, due 
to charge transfer from Fe to Rh, combined with the extremely reduced local coordination. 
The Rh local moments are important already in the pure clusters (N < 8). Therefore, they 
are not significantly enhanced by Fe doping. However, the overall stability of magnetism, as 
measured by the total energy gained by the onset of spin polarization, is found to increase 
with increasing Fe content. 
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FeRh clusters are expected to develop a variety of further interesting behaviors, which still 
remain to be explored. For instance, larger cluster should show a more complex dependence 
of the magnetic order as a function of concentration. In particular for large Rh content one 
should observe a transition from FM-like to AF-like order with increasing cluster size, in 
agreement with the AF phase found in solids for more than 50% Rh concentration. Moreover, 
the metamagnetic transition observed in bulk FeRh alloys also puts forward the possibility 
of similar interesting phenomena in nanoalloys as a function of temperature. Finally, the 
contributions of orbital magnetism and magnetic anisotropy deserve to be explored in detail 
as a function of composition and chemical order, even for the smallest sizes, particularly 
because of their implications for potential applications.— 
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